import numpy as np

# 创建系数矩阵
A=np.array([[1,2,3],[1,6,7],[1,10,6]])


# 创建常数向量
b= np.array([5,9,8])

# 解方程
x=np.linalg.solve(A,b)

print(x)